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We investigate magnetic properties of Mott-insulating phases of ultracold Bose and Fermi spinor 
gases in optical lattices. We consider in particular the F — 2 Bose gas, and the F — 3/2 and 
F = 5/2 Fermi gases. We derive effective spin Hamiltonians for one and two atoms per site and 
discuss the possibilities of manipulating the magnetic properties of the system using optical Feshbach 
resonances. We discuss low temperature quantum phases of a *^Rb gas in the F = 2 hyperfine state, 
as well as possible realizations of high spin Fermi gases with either ''Li or ^'^^Cs atoms in the F = 3/2 
state, and with ^^'^Yb atoms in the F — 5/2 state. 
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I. INTRODUCTION 

Ultracold spinor gases The seminal theory pa- 
pers by T.-L. Ho and T. Ohmi and K. Machida 01 
on spinor F = 1 Bose-Einstcin Condensates (BECs) as 
well as the experiments performed by the MIT group 
on optically trapped F = 1 Sodium condensates [3j, 
have brought a new perspective to the study of mag- 
netic systems using ultracold atomic gases. Interactions 
in bosonic systems with spin degrees of freedom host a 
wide variety of exotic phases at zero temperature and 
a dynamics clearly differentiated from the one displayed 
by scalar condensates. Recent studies involving both, 
F = 1 and F = 2 Rubidium atoms, have focused, for in- 
stance, on the rich dynamics of spinor Bose condensates 
[1, @, H- Ground-state and dynamical properties of 
F = 3 Bose condensates have also been discussed in con- 
nection to ongoing experiments with Chromium atoms 

i- 

Spinor Bose gases in optical lattices Several ex- 
perimental groups, who until now achieved spin F = 1 
and F = 2 condensates in optical traps, have already 
started, or are indeedplanning, experiments with bosonic 
spinor lattice gases [9|. Confining spinor BECs in opti- 
cal lattices offers a unique opportunity to study magnetic 
properties of matter, as a large range of tunable parame- 
ters exists, which are not accessible in solid state systems 
(for a review see [lo|). For instance, disorder, or pseudo- 
disorder, can be created in a highly controlled way in op- 



tical lattices 11 



12l |. This broadens enormously the va- 
riety of magnetic systems which can be " mimicked" with 
ultracold gases. Spinor gases also offer the possibihty 
of measuring quantum fluctuations of the total magneti- 
zation of the system employing quantum Faraday effect 
TgI . 14 1 . Another particularly interesting application of 



ultracold spinor gases in optical lattices is the engineer- 



ing of strongly correlated states, and processing of quan- 
tum information. The first steps toward these goals have 
been achieved by A. Widera et al. [l^, who investigated 
a method of measuring the coherent coUisional spin dy- 
namics in a lattice. This method has allowed for a very 
precise determination of the ^^Rb scattering lengths for 
F = 1 and F = 2 0. 

F = 1 gases in optical lattices Studies of F = 1 
systems have already been carried out by Demler's group 
17|. They have derived an approximate phase diagram 
for the case of antiferromagnetic interactions of ^'^Na. As 
in the standard Bose-Hubbard model, an _F = 1 spinor 
gas undergoes superfluid (SF) to Mott insulator (MI) 
transition as tunneling is decreased. In the antiferro- 
magnetic case in 2D and 3D, the SF phase is polar, and 
so are the Mott states with an odd number of atoms 
per site (those states are also termed as nematic) . In the 
case of even N, for small tunneling the Mott states are 
singlets, and for moderate tunneling there occurs a first 
order transition to the nematic state. 

In ID there is furthermore the possibility of a dimer- 
ized state, as in the Majumdar- Ghosh model This 
possibility was studied by Yip , who derived an effec- 
tive spin Hamiltonian for the MI state with = 1. Us- 
ing a variational ansatz interpolating between dimer and 
nematic states indicated that in a wide range of parame- 
ters the spinor ^■^Na lattice gas should have a dimerized 
ground state in ID, 2D, and 3D. This is a very interesting 
result, since so far dimer states have not been observed 
in experiments. This result has been supported by rig- 
orous studies in Ref. It was shown that under the 
same effective Hamiltonian, for a system with an even 
number of sites the ground state has total spin Stot = 0, 
while the first excited state has Stot = 2. Yip's results 
were recently confirmed by Rizzi et al. [2l|, who nu- 
merically studied the SF - MI transition in the F = 1 
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BosG-Hubbard model in ID. They found that in the low 
tunneling regime of the first MI lobe (where the effective 
spin model of Ref. works) the system is always in a 
dimerized state. Similar results were obtained by Porras 
et al. (2^ . Thus, strictly speaking, nematic order seems 
to be absent in ID in the thermodynamic limit. However, 
susceptibility to nematic ordering grows close to the bor- 
der of the ferromagnetic phases, indicating that it may 
persist in finite systems. 

Another interesting aspect, namely the possibility of 
controlling the order of the SF - MI transition by using 
appropriately polarized (lin-0-lin) laser fields to form the 
optical lattice was investigated in Refs. 2^, 24 1. Such a 
laser configuration couples the states with niF = ±1, so 
that the system becomes effectively two-component. 

F = 2 gases in optical lattices The mean field 
states of spinor F = 2 gases have been for the first time 



investigated in Refs. [25|, |26|, |27|. It is worth noticing 



that mean field states are also valid for MI states with 
one atom per lattice site, provided all atoms are described 
by the same single-particle wave function attached to a 
given site. Refs. [2^, [l^l go one step further, and apart 
from the mean field theory consider also the extreme case 
of quenched (immobile) F = 2 bosons in an optical lat- 
tice. In other words, these articles characterize possible 
on-site states for TV bosons with total spin S in the ab- 
sence of tunneling. 

After submission of the first version of this paper, Bar- 
nett. Turner, and Demler presented a beautiful and com- 
plete classification of the mean field phases for arbitrary 
F, based on 19th century method by F. Klein of solving 
quintic polynomials by the analysis of rotations of regular 
icosahedra We discuss their results in more details 
in the following. Also, very recently the effective spin 
Hamiltonians (in the first MI lobe), and quantum insu- 
lating phases of = 2 bosons have been studied by Zhou 
amd Semcnoff (29| . applying the variational principle to 
product (Guztwiller ansatz, cf. [lOl), dimer, and trimer 
states. Their results concerning efffective spin Hamilto- 
nians agree with ours. 

Spinor Fermi gases in optical lattices Obviously, 
there is an enormous interest also in Fermi gases, and in 
particular in spinor Fermi gases in optical lattices. The 
first reason is, of course, that such systems could realize 
a perfect quantum simulator of the fermionic Hubbard 
model, and thus shine some light on the problem of high 
Tc superconductivity. For spin F = 1/2 this has been 
proposed in Ref. [s^l , and recently considered with three- 
component fermions in Ref. [s^. 

Liu et al. 3J] proposed to use fermions with high F 



to realize spin-dependent Hubbard models, in which hop- 
ping parameters are spin-dependent. Such models lead to 
exotic kinds of superfluidity, such as to a phase in which 
SF and normal component coexist at zero temperature. 
W. Hofstetter and collaborators have written a series of 



symmetry in optical lattices. Such systems also have ex- 
otic superfluid and flavor-ordered ground states, and ex- 
hibit very rich behavior in the presence of disorder. 

It is, of course, inevitable to ask which atoms can be 
used to realize high-F fermonic spinor gases in optical 
lattices. The most commonly used alkali ^Li has hyper- 
fine manifolds with F = 1/2 and F = 3/2. The latter 
is, obviously, subjected to two-body losses, but as in the 
case of the F = 2 manifold of Rubidium, one can expect 
reasonably long life time in the lattice (especially in MI 
states with N = 1). Another commonly used fermion is 
a heavy alkali ""^"K, which has manifolds F = 7/2 and 
F = 9/2. These fermions are particularly useful for spin- 
dependent Hubbard models |34l |. 

There are several atoms whose lowest hyperfine mani- 
fold has F = 3/2, i.e., in those ground states two body 
losses can be avoided: ^Be, ^^^Cs, or ^'^^Ba, but so far 
only the bosonic Cesium BEC has been achieved [1^. 
On other hand, recently a BEC of ^^''Yb atoms [l^l, as 
well as degenerate gas of ^^^Yb fermions with F = 5/2, 
has been realized. Finally, fermionic Chromium has 4 
hyperfine manifolds with F = 9/2,7/2,5/2,3/2, in the 
ascending order of energies, and after achieving BEC of 
the bosonic Chromium 3l|, the prospect for achieving 



papers, reviewed in 35[, on fermionic atoms with SU{N) 



ultracold degenerate Fermi gases are very good. 

F = 3/2 and F = 5/2 Fermi gases in optical lat- 
tices Recently, there has been a lot of progress in 
understanding the special properties of F = 3/2 and 
F = 5/2 Fermi gases. C. Wu et al. HI] realized that the 
spin-3/2 fermion models with contact interactions have a 
generic 5*0(5) symmetry without any fine-tuning of pa- 
rameters, and employed this fact to propose a quantum 
Monte Carlo algorithm free of the sign-problem. They 
have found novel competing orders (39| in spin-3/2 cold 
atomic systems in one-dimensional optical traps and lat- 
tices. In particular, the quartetting phase, a four-fermion 
counterpart of Cooper pairing, exists in a large portion of 
the phase diagram. Recently, these authors studied the 
s-wave quintet Cooper pairing phase (having Stot = 2) 
in spin-3/2 atomic gases, and identified various novel fea- 
tures which do not appear in the spin- 1/2 counterpart. 
For instance, a single quantum vortex was shown to be 
energetically less stable than a pair of half-quantum vor- 
tices. 

The bosonization approach was applied to ID systems 
with F = 3/2, 5/2, ... by Lecheminant et al. [4l[. They 
used a somewhat simpler model with a spin independent 
coupling U and a coupling in the singlet channel V , and 
studied the phase diagram in the U — V plane. They 
found 3 phases: a spin density wave, an atomic density 
wave (which may crossover to a molecular superfluid), 
and a BCS superfluid (that may crossover to a molecular 
density wave). They have also classified Mott phases at 
commensurate l/{F+\/2) fiUings. Finally, the F = 3/2 
model has been solved analytically using Bethe ansatz 
and identifying the effective low energy field theory de- 
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scribing the system to be that of a deformed Gross-Neveu 
model with either 0(7) x Z2 symmetry at half-fiUing, or 
U{1) X 0(5) X Z2 symmetry otherwise [i^l- An overview 
of hidden symmetries and competing orders in spin 3/2 
gases is presented in the excehent pap er of C. Wu [43|. 
In a very recent preprint, Tu et al. |44l | have studied spin 
quadrupole ordering in spin-3/2 gases and derived effec- 
tive spin Hamiltonians for the first and second MI lobe. 
Our results agree with theirs and provide complementary 
approach and discussion. 

Goals of the paper Concerning the Bose gases, the 
main goal of the present paper is to make predictions 
for the planned experiments with ultracold F = 2 *^Rb 
atoms in the strongly correlated regime. We will thus 
concentrate on characterizing possible MI states in the 
limit of weak tunneling. In this sense we will generalize 



the results of Refs. [26|, |27|. We will disregard here the 



Zeeman effect by assuming sufficiently efficient magnetic 
shielding. Also, we will only consider the MI states with 
A'^ = 1 or iV = 2 atoms per site. Higher atom numbers 
will inevitably lead to 3-body losses, and will be thus 
much more difficult to realize experimentally. 

We will also analyze the possibility of exploring pa- 
rameters of the systems by modifying atomic scattering 
lengths. This cannot be done using the standard Fesh- 
bach resonances (cf. [i^), since we assume zero magnetic 
field. Instead, one has to use the method of optical Fes- 
hbach resonances, proposed by Fedichev et al. 
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Our main goal will be to derive effective spin models in 
the experimentally relevant regimes to study the ground 
states of the systems, and in particular to identify those 
instances when exact solutions are available, either in 
the form of product (mean field) states, or the so called 
matrix product states (MPS) l48|, similarly as it happens 
in the famous AKLT model [49|. 

Similar goals concern Fermi gases, although there 
much less is known about values of scattering lengths, 
the possibility of optical Feshbach resonances, etc. We 
will take our liberty here to explore larger regions of pa- 
rameters, assuming optimistically that they will become 
feasible experimentally at some point. 

Plan of the paper In Sec. II we discuss the Bose- 
Hubbard Hamiltonian of = 2 bosons in an optical 
lattice and its ground states when atoms are quenched 
at fixed lattice sites. In Sec. Ill we derive an effective 
spin Hamiltonian for this system. In Sec. IV we inves- 
tigate which types of ground states could be achieved, 
assuming a (limited) experimental control over the spin- 
dependent scattering lengths. In Sees. V and VI we ana- 
lyze the Fermi- Hubbard and the effective Hamiltonian for 
F = 3/2, and in Sees. VII and VIII we perform a similiar 
analysis for F = 5/2. We conclude in Sec. IX. Appendix 
A contains a detailed analysis of possibilities of optical 
manipulations of scattering lengths of *^Rb atoms in the 
F ~ 2 hypcrfine manifold. Appendix B gives a short 
overview of MPS and PEPS methods. 



II. F = 2 SPINOR BOSE-HUBBARD 
HAMILTONIAN 

The system We consider F = 2 atoms at low tem- 
peratures confined in a deep optical lattice so that the 
system is well described by a Bose-Hubbard Hamiltonian 
[30| . We assume atoms to interact via a a zero-range po- 
tential. The total angular momentum of two colliding 
identical bosons is restricted to even values due to Bose 
symmetry, so that the s-wave interaction between two 
spin-2 particles can be written asV = ffo-Po+ff2^2+54-P4, 
where Ps {S = 0, 2, 4) is the projector onto the subspace 
with total spin S. The interaction strengths, gs, depend 
on the total spin of the two colliding particles. For the 
different channels they are given by the scattering lengths 
as through gs= 4-^. 

To better understand ground state properties of the 
Hamiltonian it is convenient to express the interaction 
potential V in terms of spin operators. Making use 
of the identities / = Pq + A + A and Fi ■ F2 = 
— 6P0 — + 4P4, where Fi corresponds to the spin op- 
erator of particle i, we find V = cqI + ciFi ■ F2 + C2P0. 
Here cq = (3^4 + 4^2 )/7, ci = (54 -_52)/7, and C2 = 
(3^4 — 10^2 + 7 go)/ 7 [26|. Following [27|, Pq can also be 
expressed in terms of " singlet pair" creation and annihi- 



lation operators 5-)- = dJaJ/2 — a{aLi +020 



-2i 



s- = si 



where dj. (So-) creates (annihilates) a particle with spin 
projection a. The operator S-^- applied on the vacuum 
creates, except for normalization, two bosons in a spin 
singlet state. Such a pair does not represent a composite 
boson, as S+ and S- do not satisfy Bose commutation re- 
lations. Pq can now be written in second quantization as 
Po = 2S+S^/5. Its eigenvalues are Ns{2N -2Ns + 
where the quantum number Ns denotes the number of 
spin-singlet pairs and N the total number of bosons [i^l . 

Bose-Hubbard Hamiltonian The Bose-Hubbard 
Hamiltonian (BHH) for spin F ~ 2 can be written as 

H ^-t {al.a^j + a^a<,i) + ^ gsPsi, (1) 

<ij>,<T 

or alternatively 



i,S 



<ij>,(7 



(2) 



where 0,^1 annihilates a particle in a hyperfine state 
mp = (T at site i, Ni = X^o- '^ili^a-i is the number of 
particles at site i, Fi = ali Tc^a'aaH is the spin op- 

erator at site i ( T^^i being the usual spin matrices for 
a spin-2 particle), and : X : denotes normal ordering of 
the operator X. The coefficients Ci = Ci J d^x \w{x)\^ , 
and gs ^ gs J \w{x)\'^, where w{x) is the Wan- 
nier function centered at a; = 0. The first two terms 
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in the Hamiltonian represent tunneling between nearest- 
neighbor sites and Hubbard repulsion between atoms on 
the same site, respectively, as in the standard Bose- 
Hubbard model. The two remaining terms represent 
the energy associated with spin configurations within 
lattice sites. The ratios between the various interac- 
tions, ci/co and C2/C0, are fixed by the scattering lengths 
(ci/co = ci/co and C2/C0 = C2/C0). We assume here that 
the scattering lengths are such that the Hamiltonian ([2|) 
is stable with respect to collapse. Stability requires the 
gs > for all S, which is the case, e.g., for ^^Rb. The 
ratio t/co between tunneling and Hubbard re pul sion can 
be tuned by changing the lattice parameters [30| • When 
t <C Co, the system is in a Mott-insulating phase in which 
atoms are quenched at fixed lattice sites. Wc consider 
here the case when tunneling is sufficiently weak com- 
pared to Hubbard repulsion so that it can be treated as 
a perturbation with t/co being a small parameter. We 
study systems with one and two particles per site, which 
are the most interesting cases as they do not suffer from 
three-body losses. 

On— site Hamiltonian To zeroth order, the Hamil- 
tonian is a sum of independent single-site Hamiltonians 
(we omit the index i): 



Ho = |iV(iV-l) + | 



F F 



2C2 



S+S-. 



(3) 



Exact eigenstates of this Hamiltonian have been obtained 



in Ref. 27[. Since S± commute with the total spin op- 



erator, the energy eigenstates can be labeled with four 
quantum numbers as \N, Ns, F)mp, where is the num- 
ber of particles per site, Ns the number of spin-singlet 
pairs, and F is the total on-site spin. The eigenstates 
have a 2F + 1-fold degeneracy associated with the quan- 
tum number 1 



Their energies are: 



E = ^N{N -1) + ^{F{F+1)-6N) + 

C2 



^Nsi2N-2Ns + 3). 
5 



(4) 



In general there is an additional degeneracy which, how- 
ever, manifest itself only for states with larger number N 
of particles per site than considered in this paper. Single 



particle states are denoted as |1,0, 2), 



be- 



ing |n) the vacuum. Two and three particle states with 
maximal spin projection mp are [27| : 

12,0,2)2 



^[2V2' 

12,0,4)4 -^44lf^), 



/210 ^ 



|3, 0,0)0 = 

3\/3a\a\aL2 + 3-\/3a2aLiaLi — GV^i 
13,1,2)2 = ^ [, 
|3,0,3)s 



a2aQaL 



/20 



|3, 0,4)4 



1 

V22 

13,0,6)6 = ^4 



2\/2alalaQ 



V3alala\]\n), 



ial\n). 



States with lower mp can be obtained by acting with the 
spin lowering operator on the above listed states. 

Phases at t=0. It is easy to check which phases will 
be realized in the limit of vanishing tunneling for a given 
chemical potential. For fi < the state with no atoms has 
the smallest (Gibbs potential) energy G = E — = 0. 
For /i > we enter the phase with one atom for site 
|1, 0, 2) with G = —fi. As we increase fi further, we enter 
one of the phases with two atoms per site, namely the 
one that correspond to the smallest gg and G ^ gs ^ 2fi: 
i) |2,1,0) if .90 < 32,54, ii) |2,0,2) if 52 < ,9o,54, and iii) 
|2,0,4) if gi<g2, go- 



Ill. EFFECTIVE HAMILTONIAN (BOSONIC 
CASE) 

To derive the effective Hamiltonian to second order 
in t we consider the two-site problem. The tunneling 
Hamiltonian = —t^^ <jj> ('*L^'^i+'^ lino's) conserves 
both, the total spin S and the projection of the total spin 
[iZl- Thus, to second order, the shift of the energy of 
the two-site ground state \g, S) with total spin S is given 

by 



-E 



mHt\g,S)\^^ 



E,, — E, 



g.s 



(5) 



where v labels the (virtual) intermediate states and 
Eu, Eg s denote the unperturbed energies of the two- 
site states 1 1/), \g,S) (which are non-degenerate apart 
from the mp degeneracy). The dependence of the en- 
ergy shifts on the total spin of the two sites introduces 
nearest-neighbor spin-spin interactions in the lattice. It 
is sufficient to evaluate these shifts for only one value of 
the projection ms of the total spin. This is because tun- 
neling cannot mix states with different ms and overlaps 
Kz^ji/flg, S") I are rotationally invariant. To simplify the 
calculations we choose always the highest possible value 
of ms- 



A. One atom per site 

Pair Hamiltonian For a single particle per site, i.e., 
|1, 0, 2)'*) (g) |1, 0, 2)^^\ only 6 intermediate states are pos- 
sible: 

|f7)(*) ® |2,1,0)(^') andi^ j (5 = 0), 
|f7)« ® |2,0,2)(^') andi^j (5 = 2), 
|rj)W® |2,0,4)(^') andi^j (5 = 4). 
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The corresponding energy shifts are 



B. Two atoms per site 



4^2 
9s '' 



(6) 



or, written in terms of the c^'s, 
4t2 



it' 



eo 



Co + C2 - 6ci 



£2 



Co - 3ci 



£4 



Co + 4ci 



In this case, the overlap {i'\Ht\g,S) is the same for all 
virtual tunneling states Therefore, eg depends only 
on the difference of scattering lengths as- This suggests 
that control and engineering of the magnetic properties 
of the system could be achieved using optical Feshbach 



resonances 
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50|. One can in principle also use mag- 
netic fields to control the scattering properties, but that 
would inevitably lead to (hnear and/or quadratic) Zee- 
man effects, which would change the structure of Hq (Eq. 
([3])). Corresponding effects will be discussed elsewhere. 

The effective spin-spin Hamiltonian in second order 
reads 



H 



(^3) 



£2 A 



£4Pz 



(7) 



For ^ Rb and using the scattering lengths at zero mag- 



netic field, the energy shifts are eo 



-(4t7.9o), £2 



-0.962(4t7go), and e4 = -0.906(4tV5o). Thus, as eo is 
smallest, ^^Rb should experience antiferromagnetic be- 
havior in the ground state (cf. [l^). For ^'^Rb, the scat- 
tering lengths are negative in the absence of a magnetic 
field. Thus Hamiltonian ([2]) is unstable with respect to 
collapse. One can, however, use optical Feshbach reso- 
nances (similarly as has been demonstrated using mag- 



netic Feshbach resonances l51|) to achieve gs > for all 
S. In this case we expect the physics of the ®^Rb lattice 
spinor gas to be similar to the case of ^^Rb. Despite the 
fact that control of magnetic properties is possible using 
optical Feshbach resonances, observation of the ground 
states requires, in turn, very low temperatures to resolve 
accurately the different energy shifts. 

Lattice Hamiltonian The Hamiltonian ([7|) can be 
easily generalized to the whole lattice, H = J2i ^o.i + 



E 



<ij> 



(y) 



It can also be transformed into a polyno- 



mial of fourth order in the Heisenberg interaction • Fj 



39eo - 80e2 
51 



(F. • F,) 



(F.-F,f 



9eo - 8e2 
102 

7eo 10e2 64 
204 204 72 
7eo + 10e4 



(F, • Fj)' 



1020 



-(F.-F,r 



(8) 



Here up to four powers of Fi ■ Fj appear due to the fact 
that channels 1 and 3 between different sites are not for- 
bidden. 



Accessible single site states We now derive the 
effective Hamiltonian for the case in which the ground 
state corresponds to two particles per site and proceed 
as before by considering the two-site problem. To cal- 
culate the energy shifts, we find first the ground states 
of the unperturbed Hamiltonian Hq, which are easily de- 
termined from Eq. (|4]). The ground state corresponds 
to 

• |2, 1,0) if go < 32,54, «-e., for ci < 0, C2 < lOci and 
ci > 0, C2 < 3ci. This singlet state is the single- 
site ground state of ^^Rb for unmodified values of 
scattering lengths. 

• |2, 0, 2) if g2 < go, g4, i-e., for ci > 0, C2 > 3ci. This 
situation may be accessed with optical Feshbach 
resonances (see Appendix A for details). 

• |2, 0, 4) if 54 < g2, go, *-e-, for ci < and C2 > lOci; 
this case is hardly accessible experimentally. Using 
optical Feshbach resonance to achieve it, would lead 
to enormous losses (see Appendix A). 

In the case of singlets |2,1,0), the lattice ground 
state is non-degenerate in zeroth order in t, such that 
it does not have any effective dynamics. It does, how- 
ever, have the first order correction to the wave func- 
tion I*) = (1 - {Ho - i;o)-i^t)nj2,l,0). It has 
also the second order shift of the ground state energy 
SE^-(^P\Ht{HQ-Eo)Ht\ij). 

Pair Hamiltonian To calculate the energy shifts to 
second order, we assume for simplicity that the energy 
spacing between eigenstates of Ho is sufficiently large 
compared to tunneling transitions to intermediate states 
and treat; therefore, these states as non-degenerate. We 
shall consider here only the |2,0,2) single-site ground 
state, because the state |2, 0,4) is hardly accessible ex- 
perimentally. The corresponding effective Hamiltonian 
admits a very large number of virtual intermediate states 
and is, therefore, very complex to calculate. We ex- 
pect, however, that the physics for the cases |2,0, 2) and 
1 2, 0, 4) is similar. 

Starting from |2,0,2)W (g, |2,0,2)(-'\ to second order 
in t there are 26 intermediate virtual states spanning the 
five channels of total spin S* = — 4: 

|1,0,2 >W (g)|3, 1,2 >(^^ and i '-^ j (5 = 0.. .4), 

|1,0,2 >W ®|3,0,0 >(^') and i <^ j (5 = 2=, 

|1,0,2 >W ®|3,0,3 >(^^ and i j (5 = 1.. .4), 

|1,0,2 >W (g)|3,0,4 >(^') and i '-^ j (5 = 2.. .4). 
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The corresponding energy shifts are: 



eo 



ei 



1 (—f 

4.2 ^ V 35 ; 



7 Co + 7/5c2 ' 

1 (ir 



2 m" 



7 Co + 7/5c2 21 Co + 3ci 



£2 = -t^ 



1 



1 (W)' , 3(1)^ 



1 



7co + 7/5c2 15 Co — 3ci cq + 3ci 



' 1124 ^ 
. 49 ; 



11 • 35 Co + 7ci 



£3 = 



£4 



13 



7co + 7/5c2 7co + 3ci 11 co + 7ci 



7co + 7/5c2 Co + 3ci cq + 7ci 



The full effective spin Hamiltonian reads H = i/o i + 

E<.,>(^oPo^^^VeiA^'-'V£2Pi^^Ue3Pi'^Ue4Pf ^ If 
scattering lengths are changed through optical Feshbach 
resonances to have |2,0,2) as the unperturbed single- 
site ground state, H has a ferromagnetic lattice ground 
state l^-) = 0. |2,0,2)^'^ for ^^Rb. The energy shifts of 
the different spin channels arc now mostly determined 
by the different contributions of the intermediate states 
\v) and Clebsch-Gordon coefficients rather than by the 
differences in scattering length. As a result we do not 
expect that optical Feshbach resonances will allow such 
a precise and extensive control of spin-spin interactions. 
We expect that in the |2, 0, 4) phase energy shifts exhibit 
a similar behavior determined essentially by the values 
of Clebsch-Gordon coefficients, although the magnetic 
properties of this phase have yet to be calculated. Quite 
generally, we conjecture that tuning spin-spin interac- 
tions via optical Feshbach resonances in systems with two 
particles per site will be more difficult than for systems 
with one particle per site. 



IV. GROUND STATE PROPERTIES 

In this section wc will study the ground state proper- 
ties of systems with one or two particles per site, con- 
centrating on the question whether the control over the 
scattering lengths can lead to the appearance of some 
more exotic phases (c.f. 53), assuming a simple one di- 
mensional chain or a two-dimensional square lattice. To 



this aim we use the approach developed recently by Wolf 
et al. [Hst . and search for such combinations of parame- 
ters for which we can represent the ground state of our 
Hamiltonian exactly (isolated exact ground states), or 
nearly exactly using matrix product states (MPS) [48| in 
ID, and projected entangled-pair states (PEPS) in 2D 

M. 



First we add to the bond Hamiltonian ([7]) a certain 
number of times the identity operator J^*-'^ = J2s ^s'''^ 
on the bond, so that the Hamiltonian becomes positive 
definite, i.e., 



H 



s=o 



(9) 



with all As being non-negative. 

In the case of one atom per site (see Subsection A) 
Al = A3 > take the greatest values, whereas one of the 
Xs for = 0, 2, 4 can be set to zero. We will subsequently 
assume that the control over the scattering lengths per- 
mits to choose which one is set to zero and to fix the 
magnitude of the others, even though this goes beyond 
the experimental feasibility; we will treat experimentally 
accessible cases with particular attention. 

In the case of two atoms per site, the values of As 
are essentially determined by the Clebsch-Gordan coeffi- 
cients. They are descending quite significantly as S in- 
creases, with A4 = 0. Thus an extensive control of the 
effective Hamiltonian using optical Feshbach resonances 
is hardly possible. Nevertheless, in Subsection B we ex- 
plore some limiting cases by setting further As's to zero. 

Seaching for isolated exact ground states Our 
approach to search for exact (or, at least "variation- 
ally exact") ground states can be described as follows. 
In ID with periodic boundary conditions we seek for 
translationally invariant MPS of N spins F = 2, given 

by I*) = EH,..,...,^„Tr(A^^^'^...A'«)|zi,Z2,...,*^), 
where i's enumerate the computational spin 2 basis, 
while A's are d x d matrices, d > 1. Let R{ij) de- 
note the range of the reduced density matrix pij = 
Trfc^,j(|«')(«'|). This range is given by R{ij) = 
spaujv/ Tr(M^M^)|z, j), where the M's are arbi- 
trary d X d matrices. Let K{ij) denote the kernel of 
the bond Hamiltonian Note that if R{ij) C K{ij), 
then H\'^) — 0, and since H is non-negative, that im- 
plies that I^P) is a ground state. Similarly, in 2D with 
periodic boundary conditions we seek for translation- 
ally invariant PEPS of spins = 2 in 2D; given 
by I*) = E.„.,,...,.„ Cont(A»M'^ . . . A^")|zi, Z2, . .' . , zat), 
where ik's enumerate the computational spin 2 basis of 
the k-th atom, and k denotes coordinates in the 2D 
square lattice. This time the tensors A^'' with four 
indices, all running from 1 to d, correspond to the 4 
bonds touching the fcth site. The contraction is over 
pairs of indices belonging to the same nearest neigh- 
bor bond. In this case R{ij) is given by R{ij) — 
spauj^j^- ■Cont(M A'- A^)\i,j), where the M's arc ar- 
bitrary tensors with the 6 indices that are not con- 
tracted in the product of A^A^ . Alternatively, we may 
search for states that break the translational symme- 
try; antiferromagnetic Neel-like ordering in ID could 
for instance correspond to MPS of the form Y^) = 

T.^,,^,,...,^^^M''B'-A'^ ...B'-)\i^,12,...,In) for CVCU 
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N. 



A. One atom per site 

Mean field diagram Before we proceed, it is worth 
to discuss the mean field phase diagram obtained un- 
der the assumption that the ground state is a product 
state, l^*) = |e, e, . . .) (see [25|,l28|). We are following here 
the most complete description of the phase diagram, pro- 
vided recently in (28| . We introduce the nematic tensor 

Qab - ^{Fah + FbFa) " ^SabiF^). Thcrc arc 3 possible 
mean field {i.e., product) ground states, with |e) given 
up to 5*0(3) rotations: 

• Ferromagnetic state, |e) = (1,0,0,0,0); possesses 
only the U{1) symmetry of rotations around the z~ 
axis, and has maximal projection of the spin onto 
z axis. 

• nematic state, which apart from 5*0(3) rota- 
tion, has an additional 77-degeneracy, |e) ~ 
(sin(77)/\/2, 0, cos(r/), 0, sin(?7)/v^). This state is a 
MI version of the polar state in BEC; it has mean 
value of all components of the spin equal zero, but 
non vanishing singlet projection (singletje, e) 7^ 0; 

• Tetrahedratic (cyclic) state, |e) = 
(1/a/3,0,0, \/273,0); this is a MI version of 
the cyclic state. The state may be uni- or biaxial, 
depending on whether the nematic tensor does, or 
does not have a pair of degenerated eigenvalues; it 
has vanishing of both, of mean values of all of the 
spin components, and of the singlet projection. 

The phase diagram is such that the system is in: 

• a ferromagnetic state for A4 = 0, A2,Ao > 0, and 
for Aq = 0, provided A2 > I7A4/IO; 

• a nematic state for Aq = 0, provided 3A4/IOA2 < 
I7A4/IO; 

• a cyclic state for Aq = 0, provided A2 < 3A4/IO, 
and for A2 = 0. 

We will identify below the regimes of the phase dia- 
gram in which the mean field diagram is exact, or nearly 
exact. 

MPS reduce to mean field states Since Ai ~ 
A3 > take the greatest values, the ground state of ([9|) 
should belong to the symmetric subspace. For transla- 
tionally invariant MPS this implies [A*,^-'] = 0, i.e., 
the same matrix K transforms A* and into the 
commuting Jordan forms. Generically, if all eigenval- 
ues A^, A-[ of A^,A^ are distinct, such MPS correspond 
to linear combinations of product states |efc,efc,...), 
where |efc) = X^i^^li). Since the Hamiltonian is a 
sum of nearest neighbor bond Hamiltonians, we have 



I]fc,fc'('Sfc,efc . . . |iy7|efc/,efc' . . .) cx (efcle^/)^ vn z. ID- 
chain, and thus in the limit of an infinite chain the 
ground states are equally well described by product states 
(that will typically break the rotational symmetry) . This 
means in this case we expect mean field (product) states 
to provide a very good approximation of the ground 
states with translational symmetry. Let us analyze the 
different possible cases: 

(Ai) For A4 = A2 = Ao = 0, all symmetric states 
are ground states, i.e., in particular all product states 
|e, e . . .) with arbitrary |e). 

(Bi) For A4 = A2 = 0, Aq > 0, the ground states 
|e, e . . .) remind the cyclic states states of Ref. [13] [i.e., 
they correspond to translationally but nor rotationally 
invariant product states), which now mix 5 = 2 and 
5* = 4 contributions on each bond, and they have to 
fulfill the condition (singletje, e) = 0. Denoting by |e) = 
(e2, ei, eo, e_i, e_2), this implies Bq— 2eie_i-f 2e2e_2 = 0. 
These states form a much greater class than the cyclic 
ones, since they may have non- vanishing (and even maxi- 
mal) components of the spin. Interestingly, the transition 
between the cyclic phase for A2 = 0, and the ferromag- 
netic phase for A4 = 0, occurs via such states, i.e., at 
the transition point the degeneracy of the ground states 
manifold explodes. 

(Ci) For A4 = and A2, Aq > 0, the ground states are 
ferromagnetic states |2)n|2)n ■ ■ • |2)„, corresponding to a 
maximal projection of the local spin onto a given direc- 
tion n= (sin(0) cos((?i)), sin(0) sin((/)), cos(0)). Such vec- 
tors for F = 2 may be parametrized (in the basis of Fn 
with descending rap) as cx (z~^, 22;"^, -\/6, 2z, z^) with 
z = \z\e^'^ , \z\ € (—00,00). It should be stressed that 
ferromagnetic states are exact ground state in the entire 
part of the phase diagram whenever A4 = 0. 
(Di) For Ao = and A4, A2 > 0, the ground states 
apparently favor antiferromagnetic order. This, how- 
ever, can be misleading, if A4 ^ A2. In that case, as 
the mean field diagram suggests, the ferromagnetic or- 
der might prevail. We have applied in ID a more gen- 
eral variational approach, going beyond mean field. We 
have looked for ground states by applying the variational 
principle to mean field (product) states |e, e...), Ncel- 
type states |e, /, e, / . . .), and valence bond solid states 
with singlet states for distinct pairs (dimers) of neigh- 
boring atoms and translational dimer symmetry. For 
the mean field case as discussed earlier the energy is 
either minimized by the ferromagnetic state |e) = |2)n 
(for A2 > I7A4/IO), by a nematic state |e) = |0)n (for 
3A4/IO < A2 < I7A4/IO; in this case the state is a com- 
bination of total spin 0, 2 and 4), or, for A2 < 3A4/IO, 
by a cyclic state, |e) = (e2, ei, eo, e_i, 6-2) with eo = 
l/-\/2, e2 = — e_2 = 1/2, e\ = e_i = 0. Imposing Neel 
order with (e|/) 7^ 1 always results in a larger energy, 
as Ai,3 > A2,4, and the overlap with the singlet can be 
maximized already by restricting to product states. On 
the other hand, for the dimer state the energy per bond 
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Figure 1; Sketch of the phase diagram for the (Di) case, 
obtained by applying the variational principle in the A2, A4 
phase space (for Ao ~ 0) to mean field, Neel, and dimer 
states with one atom per site. The scale is set by letting 
Ai = A3 = 1. Neel- type states are never favorable over ne- 
matic states. The ferromagnetic region (gray) was obtained 
numerically by imaginary time evolution of MPS, and com- 
paring the result from runs with d = 1 and d = 5 in a chain 
of 50 sites with open boundary conditions [s^. Of course, on 
the line (A4 = 0, A2) ferromagnetic states give always ground 
states. Dashed lines indicate the regions where the type of 
mean field state with lowest energy changes qualitatively (see 
text for more details). The red (dashed-dotted) line gives 
the values of (A2, A4) which can be obtained by changing the 
spin-independent scattering length cq — (834 +4g2)/7 of *'^Rb 
through optical Feshbach resonances. The arrow gives the 
values for unchanged co, black and white circles indicate a 
change of Co of 10% and 100%, respectively. 



is given by iTr(7J/^ 



25 - 



1). This results in the phase 
diagram shown in Fig.^^ somewhat analogous to the re- 
sults obtained by Yip [l9| . The red line depicts values of 
A2, A4 experimentally accessible through modifications of 
scattering lengths. We have applied MPS code to search 
numerically for the exact ground states using the method 
of [5^. We confirmed that in the shaded region in Fig.[T] 
there is indeed a ferromagnetic ground state. We have 
also studied the ground state at the experimentally ac- 
cessible line, and found indications of nematic and dimer 
order in the phase diagram, see Fig.[T] We expect that in 
2D, in addition a possible ground state could be formed 
from resonating valence bond states [H^l , and we are plan- 
ning to apply the 2D PEPS methods to investigate this 
question. 

(El) For A2 = and A4, Aq > 0, as in the (Di) case, 
mean field cyclic states are favorable over Neel states. 
We have compared them variationally to the analogues 
of the dimer states in the present case, i.e., configura- 
tions which have a state with total spin 5 = 2 on dis- 
tinct bonds. We call these state para-dimers. Now the 




Figure 2: Analysis of the ground states obtained variation- 
ally from an MPS ansatz for the case (Di), for combinations 
of A2 and A4 lying on the dotted-dashed line of Fig. [1] (a) 
Singlet projection SP = Xlitr(PoPi ® Pi), where pi is the 
reduced density matrix of site i, (b) nematic order param- 
eter Q = max„ Qn = maxn 5]lJ(nnSi)^ — 2]/N (nn is the 
unit vector pointing in direction SI), and (c) dimerization 
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AJV/2,iV/2+l 




pN/2- 



l,N/2 



) . (Values obtained for open 



chains of 16 sites, with d up to 30. 



situation is quite different from the dimerized states dis- 
cussed in (Di), as the states on the bond are not unique 
and states with different total spin projection Ms=2 can 
form superpositions. In the subspace of states having a 
para-dimer on each second bond, the Hamiltonian can 
be written as an effective interaction between neighbor- 
ing para-dimers, {H'^ii)(ij) = Yl's'=o^S'{P')^s''^ , where 
i, j now enumerate the para-dimers, and {P')s' projects 
onto the subspace of two para-dimers with total spin S' . 
For A2 = 0, Ao, A4 > 0, always X'2 < X'q, X'4, and the 
optimal superposition of para-dimers with different pro- 
jections is again the cyclic combination. On the other 
hand, given that is the lowest coefficient, combining 
neighboring para-dimers to states with total spin S" = 2 
might lead to even lower energies. On this level again 
states with different total z-projection Ms' can be com- 
bined, and it turns out that again the cyclic combination 
minimizes the energy. Comparing the energies of cyclic 
product states, cyclic states of para-dimers, and cyclic 
combinations of "para-dimerized" para-dimers, the phase 
diagram shown in Fig. [3] is obtained. 



B. Two atoms per site 

We discuss here cases where the values of the As 
are in descending order, determined essentially by the 
Clebsch-Gordan coefficients. In such a situation only a 
ferromagnetic-like order in the ground state is possible, 
but now the states do not neccessarily have to be of the 
product form, especially if Ai 0, or Ai = A3 = 0. We 
explore the following limiting cases: 
(A2) For all A5 = except Aq > 0, the ground states. 
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Figure 3: Sketch of the phase diagram for the case (Ei) 
case, i.e., when A2 is the lowest coefficient. As in Fig. [T] 
the Ifamiltonian is shifted and rescaled, such that A2 = 0, 
< Ao, A4 < i, and Ai = A3 = f . The lines give the bound- 
aries obtained from variationally comparing product states 
(where cyclic states are always optimal) to product states of 
para-dimers (where again the cyclic combination gives min- 
imal energy), and combinations of neighboring para-dimers 
with total spin 5" = 2. The red (dotted-dashed) line indi- 
cates the combinations of (A4,Ao) which can be obtained by 
changing co through an optical Feshbach resonance. Black 
and white circles indicate a change of co of —10% and —50%, 
respectively. 



when reduced to neighboring sites, are either of the type 
|2)n|e) (or |e)|2)„) with |e) = (e2, ei, eo, e_i, 0) in the S„ 
basis) or in the form |l)n|e) with |e) = (e2, ei, cq, 0, 0) 
(or|g)|l)„). 

(B2) Similarly, for A2 = A3 = A4 = and Ao,Ai > 0, 
ground states correspond to product states which, when 
reduced to neighboring sites, have either the form |2)n|e) 
(or |e)|2)n) with |e) = (e2, ei, eg, 0, 0) written in the Sn 
basis, or the form |l)n|l)n- 

(C2) For A3 = A4 = and all other A5 > 0, the ground 
states are product states formed by the vectors |2)n, and 
I l)n7 with the constraint that there arc not two 1 1)„ states 
in the neighboring sites. 

(D2) for A4 = and all other A5 > 0, the ground states 
arc as in the case (Ci) before, i.e., they are of the form 
|2)„|2)„---|2)„. 



V. F = 3/2 FERMI-HUBBARD HAMILTONIAN 

The system Let us now turn to the discussion of 
the spin-3/2 (in this and the next section) and -5/2 (in 
Sees. VII, VIII) Fermi lattice gases. The total wavefunc- 
tion of the fermions has to be anti-symmetric, implying 
that the spin of two colliding fermions can only be even. 



Interaction for two fermions with spin F in the s-wave 
channel can be written in the form 

y = .9oA) + .92A + ■•• + giF-lPlF-l. (10) 

where Pg is the projection operator on the subspace with 
total spin S and is the interaction strength, which 
depends on the scattering length (as) gs = ^'^^"^ ■ 

On-site Hamiltonian For two spin-3/2 particles 
with anti-symmetric spin-wavefunction, we can use the 
identities I = Pq + P2 and Fi • F2 = 7o-Fb + 72A to 
express the interaction in the form V = CqI + C2F1 • F2. 
Here, is the spin operator of the particle i and 7„ = 
[n{n + 1) - 2F{F + l)]/2. The total Hamiltonian in the 
limit of vanishing tunneling (t = 0) for i<" = 3/2 is a sum 
of single-site Hamiltonians of the form (omitting site in- 
dices) 

Ho =y ^ al^a^ija0aa+ 

C2 V- t t 

0.(3^(7 

where cq = (-go + ^92)/^ and C2 = (-50 + 52)/3, and 
symbols without bars are related to those with bars in the 
same way as in the bosonic sections. In the summation, 
greek letters are spin indices. 

The eigenvalues of the Hamiltonian Hq read 

E\N, F) = icoiV(iV - 1) + ic2(F(F -f 1) - ^7V), 

where N is the number of particles per site and F is the 
total on-site spin. As there are four accessible states per 
site, corresponding to the spin projections, the maximal 
number of particles is = 4. The energies for different 
values of A^ and F are listed below: 



E\N ^ 


1,^ = 


3/2) 


= 0, 


E\N = 


2,F = 


2) = 


3C2 

Co - — = .92, 


E\N ^ 


2,F = 


0) = 


15c, 
Co ^ = .90, 


E\N = 


3,F = 


3/2) 


= 3co ^ = {go + 532)72, 


E^{N = 


4:,F = 


0) = 


15C2 

6co = .90 + 5.92- 



(12) 



Phases at < = The actual ground state (GS) in the 
case of vanishing tunneling is determined by comparing 
the (Gibbs potential) energies G = Eo{N,F) - fiN of 
the above listed states. The resulting phase diagram is 
already quite complex and depends on the values of co 
and C2, or better to say go, 52, in a nontrivial way (see 
Fig. [5]). Denoting the states by {N,F), and writing S = 
92/50 and p, = ^/go, we obtain 
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• (0,0) is the GS for fl< 0; 

• (1,3/2) is the GS for < /!< 1, and /2 < 5; 

• (2,0) is the GS for 1 < ^ < (5(5 - l)/2, and 5 > 1; 

• (2,2) is the GS for (5 < 1, and (5 < /l < {35 + l)/2; 

• (3,3/2) is the GS for max[(3(5 + l)/2, (5(5 - l)/2] < 
fi < (5(5 + l)/2; 



(4,0) is the GS for (5(5 + l)/2 < fl. 
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Figure 4: Ground states for f = 3/2 for the case of no tunnel- 
ing, in the space of (5 = g^jga and /7 = m/So- The phases are 
labeled by (A'^, F) , where A'^ denotes the number of particles 
per site and F their total spin. 



VI. EFFECTIVE HAMILTONIAN FOR F = 3/2 

We foUow here the same lines as in bosonic part. 
We assume weak tunneling, such that we can use 
perturbation theory to calculate an effective spin- 
spin Hamiltonian. The tunnehng Hamiltonian Ht = 
— tX]<y> (t('^L%o- + O'la'^io-) before involves nearest 
neighbors only and conserves the total spin S as well as 
its z-component Ms- We apply Eq. ([5]) to calculate en- 
ergy shifts to second order in t, limiting the calculations 
to one value of Ms, which we will choose to be the largest. 



A. One atom per site 

We consider first the case of one particle in each lattice 
site. For t = 0, the ground state of the nearest neighbour 
pair with total spin S and maximal total ^-projection can 
be written as 

\SMs = 5),,, = J2 {FFm,m,\SMs = S) 4rnAnJ^)^ 

7711 ■'>^^'y 

(13) 

where is a spin of a fermion (mi, m2 are the z- 
components of this spin), and i and j are lattice indices. 



Possible intermediate states are those having two par- 
ticles on one site (say i), and no particles on the other 
(say j). Two-particle states with total on-site spin F and 
maximal projection Mp = F read 



= \F,Mp = F),, 

{FF 1711171211^, mp — F) a. 



V2 



m) 



From anti-commutation relations for fermions and prop- 
erties of Clebsch-Gordon coefficients, it follows that F 
has to be even. The energy shifts thus are simply 



4*2 

9s ' 



es= 0, 



(total spin) S is even, 
S is odd. 



(15) 



This in fact is independent of the spin F. The explicit 
expressions for F = 3/2 can be written 



£2 



Co - 3c2/4 
4*2 



4*2 

92 

4*2 

go co-15c2/4 

£3 =£1 = 0, 



(16) 



and the effective Hamiltonian is Hj'''" = £2^*2 "' + £4^4 
This result is the same as discussed in Ref. The 
effective spin model has several particularly interesting 
limits. In particular, when gg — g2 the model has a 
5 1/ (4) symmetry, and is integrable via Bethe ansatz in 
ID [45[; its ground state is a spin singlet with gapless 
excitations. 



B. Two atoms per site 

We now consider the case of two spin- 3/ 2 fermions on 
each lattice site, so that the total spin on-site may take 
values or 2. We express two-site states in the form 
\S,F), where S is the total spin of the two sites and F is 
the total on-site spin. Again we can limit calculations to 
the maximal total z-component Ms = S. We write the 
ground state in the same form as in Eq. ()13p 



\S,F), 



- {FFMN\SMs = S){FFmim2\FM) ■ 

rai ,7712 
ni ,712 , -A/, TV 

• {FF n,n2\FN) a^, aL. «Lh 1^^)^ (17) 



where M , N are the z-components of the total spins of 
sites i, k, respectively. Intermediate states \v) now have 
three fermions in site i (which is equivalent to having a 
spin-3/2 hole) and one in the neighboring site k. As in 
Eq. P?|) . we explicitly give the intermediate states with 
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total spin S and maximal total z-component Ms- Such 
virtual states exist only for integer S* = . . . 3. 



|5 = 3) 
\S = 2) 

\S=1) 
\S = 0) 



\n) 



2a^3 a^i 



%/3aJ 3 a|_ 1 3 3 ) I ) 



1 / t t t t 



a^3 , 



(18) 



Using those states, the energy shifts to second order can 
be calculated. If the total on-site spin is = 0, i.e., the 
on-site states are singlets, then those do not interact in 
second order perturbation theory, but there is a second 
order shift which amounts to eo = — 2i^/(co-l-15c2/4) per 
bond. If the atoms on-site form a composite with spin 2, 
then the energy shifts are as follows: 



£4 = £2 =0, 
£3 = £1 = - 

eo = 



4f2 



C0-9C2/4' 

10^^ 
Co - 9c2/4' 



(19) 



The above results agree with those obtained recently in 
Ref. [H. 



C. Three atoms per site 

For three particles per lattice site, possible two-site 
states are similar to the case of one atom per site, 
Eq. because states of three atoms per site can equiv- 
alcntly be written as a single hole in the filled Fermi sea. 
Intermediate states now have four particles on one site. 
Because of Pauli's principle there is only one such state, 
namely the filled Fermi sea 3/2*^! i/2'^I -1/2*^1 -3/2l^)' 
It is thus clear that the energy shifts have to be as in 
Eq. dUl): 



£2 

£0 
£3 



_ 4^2 

92 

_ ^ 

90 

:£l =0. 



4t' 



Co 



3C2/4' 
- 15C2/4' 



(20) 



Obviously, for four particles per site the on-site ground 
state is the filled Fermi sea. This is an exact eigenstate 
of the full Hamiltonian, as no tunnelings are possible in 
this state. 



VII. F = 5/2 FERMI-HUBBARD HAMILTONIAN 

The theory for insulating states of a spin-5 /2 gas with 
1 or 2 atoms per lattice site is essentially the same as in 
the case of spin-3/2 particles, so we will comment only 
the basic differences. 

The system and on-site states Now the two- 
particle s-wavc interaction can be written in the form 

V ^ goPo + g2p2 + giPi- (21) 

or, using the identity operator / and spin operators F^, 

y = coi + ci(F, •F,) + C2Po, (22) 

where cq = (5^2 + 2354)/28, Ci = (-52 + 54)/7, and 
C2 = (7^0 ~ IO52 + 334)/?. The singlet projection oper- 
ator Pq can be represented via creation and annihilation 
operators as 

Po = A^A, 

A — = ( 05 a_ 5 

\/3 V 2 2 



(23) 

a3a_3 + aia_i ] . (24) 



V3 

The on-site Hamiltonian attains then a similar form as 
in the case of spin-3/2 with an additional term C2P0 (and 
relations between Ci and Ci as before). The energies for 
different numbers of fermions per site A'^ and total on-site 
spin F arc listed below: 
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(25) 
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Phases at t = The actual ground states are deter- 
mined by minimizing the Gibbs energy G = Eo{N, F) — 
fiN. The resulting phase diagram is three-dimensional 
and quite complex, as it depends on the values of cq Ci, 
and C2 (or better to say go, gi, and 52) in a highly non- 
trivial way. The ground states are plotted in the space 
oi ^/cq and ci/cq for two values of C2 in Fig. O 

VIII. EFFECTIVE HAMILTONIAN FOR = 5/2 

A. One atom per site 

The effective Hamiltonian to second order in the tun- 
neling amplitude t has the form Hj*"'' = ^oPq''^ + 
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Figure 5: Ground states for F = 5/2 for the case of no tun- 
neling, in the space of ci/co and = ^i/(?o for cijc^ = 1/2 
(top) and C2/C0 = — 1/2 (bottom), taking into account states 
of up to 3 particles. The phases are labeled by (A'', where 
TV denotes the number of particles per site and F their total 
spin. 
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Figure 6: (a) Classification of the effective Hamiltonians 
which can be obtained for _F = 5/2 and a single particle per 
site (see Eq. in the (ci/co, C2/C0) space (for /.i/cq = 0.2). 

Thick lines indicate the borders of the region where the * = 
ground state has a single particle per site, thin lines indicate 
the borders between different regimes of the corresponding ef- 
fective Hamiltonian. When £4 is lowest, then the ground state 
is ferromagnetic. The area with eo < £2, £4, is shown in more 
detail in (b). The different regions are obtained from a varia- 
tional ansatz similar to Sec. IV, case (Di), comparing dimer 
with different mean field type states. Lowest energy states 
according to this ansatz are mostly dimerized (black). The 
ferromagnetic (dark grey), nematic (grey), and cyclic (light 
grey) regions are much smaller. 



be written in the form 



As usually 64 is smallest, the ground states in this case 
are mostly ferromagnetic, as can be seen from Fig. [5] (a). 
There is however a region where eo < £2, £4- In this case 
the variational approach followed in Sec. IV, case (Di) 
shows that within this region again dimerized as well as 
ferromagnetic, nematic, or cyclic phases can be realized 
(see Fig.EKb)). 



B. Two atoms per site 

On-site ground states with two fermions per site can 
have total spin 0, 2, or 4. Tunneling carries over these 
states into those with three atoms on one site, and one 
on the neighboring site. The state of three particles on 
lattice site i with total spin F and z-projection M can 
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As the virtual state necessarily has to be anti-symmetric, 
there are three different possibilities for F, namely 9/2, 
5/2, or 3/2. Let us now consider ground states with 
different on-site spin F separately. The case F = leads 
only to a ground state energy shift, but not to interesting 
dynamics in second order perturbation theory. 

On-site ground states with spin 4 In this case 
nearest-neighbor pairs can form total spin from S ^ 

to 5 = 8. The effective Hamiltonian in second order 

(ii) 

perturbation theory can be written in the form = 
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Tuning ci/cg and C2/C0, typically cither eg and egi £7j 
or £4 arc the smallest coefBcients, such that, though of- 
ten ferromagnetic ground states are realized, also models 
preferring anti-ferromagnetic order are possible. 

On-site ground states with spin 2 Finally, we 
consider the case of on-site ground states with spin 2, 
where 
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low temperatures. We have discussed two experimen- 
tally relevant cases with either one or two atoms per 
lattice site. Our analysis shows that in the case of a 
single atom per lattice site, the spin-spin couplings and, 
therefore, the magnetic properties of the system can be 
precisely manipulated using optical Feshbach resonances. 
We have explored the quantum phase diagram for such 
a case using variational and numerical techniques. On 
the other hand, the manipulation of the magnetic prop- 
erties of a ground state with two atoms per lattice site 
becomes much harder to achieve. In this last case, the 
spin-spin interactions present in the effective Hamilto- 
nian couple many different virtual states. Since (at zero 
magnetic field) the scattering lengths as are very similar, 
the couplings depend very strongly on the correspond- 
ing Clebsch-Gordan coefficients. In this respect, spinor 
condensates with non-alkaline atoms 3l| which present 
large differences of their scattering lengths as could dis- 
play stronger spin-spin interactions effects [81]. As it was 
pointed out to us by L. Santos, in the limit when \t/go\ 
is very small, the effective spin-spin interactions might 
involve magnetic dipole-dipole interactions. Still, for Ru- 
bidium the effective model here is valid for a certain range 
of \t/go\ for which the magnetic dipole moment can be 
neglected. On the other hand, tuning the system into 
a range where dipole-dipole interactions are important, 
they might offer an additional knob to control the effec- 
tive Hamiltonian. 

In the second part we have performed a similar analy- 
sis of F = 3/2 and F = 5/2 gases. Also there, while the 
phases with one fermion (or one fermionic hole, i.e., 2F 
fermions) per lattice site can be more easily controlled 
with optical Feshbach resonances, the physics of phases 
with 2 or 4 atoms is controlled by the Clebsch-Gordon co- 
efficients. The latter situation might still lead to various 
spin models preferring either ferromagnetically or anti- 
ferromagnetically ordered ground states. In this context 
it would be particularly appealing to realize (e.g., by 
the above mentioned dipole-dipole interactions) Hamil- 
tonians with dominant contribution of Pa-terms. Such 
Hamiltonians admit AKLT-like gapped ground states in 
2D in the honeycomb lattice [t^ . 



Now the effective Hamiltonian is i/cff = J2i ^o.i + 
E<^J> eoPo' + eiPl' + e2p;' + esP^= + e,P^ . Typically, 
either ei or eo are the smallest coefficients, such that real- 
izable spin models usually have ground states preferring 
anti-ferromagnetic order over ferromagnetic one. 



IX. CONCLUSIONS 

Summarizing, in the first part of the article we have 
analyzed the different Mott insulating phases of repul- 
sive spinor F — 2 bosons confined in optical lattices at 
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APPENDIX: OPTICAL FESHBACH 
RESONANCE FOR F = 2 RB ATOMS 

Optical modifications of scattering length, or in other 
words o ptic al Feshbach resonances (OFR) were proposed 
in Ref . [47[ , and carefully analyzed theoretically in a se- 
ries of papers by J.L. Bohn and P.S. Julienne [5^. These 
authors have pointed out that OFRs are inevitably asso- 
ciated with spontaneous emission losses, since molecular 
states used for OFR cannot be to far from the photoas- 
sociation resonance. For these reason, changes of natural 
scattering length of ^^Rb, which itself is of order 100 a.u., 
by more than 10 a.u. were considered to be unrealistic. 
These prediction have been confirmed in the recent ex- 
periment of R. Grimms group [HO] ■ 

For the present investigations this implies that only 
limited changes of scattering length are possible. Note 
that since OFR takes place far from the nucleus, where 
the excited state potential has a dipol form, —C^/R^, 
One can only modify in this way a spin independent part 
of the scattering, T hat is (804 -I- 4a2)/7. That means 
that ao and 4a4 — 3a2 remain unchanged under OFR. 

The most accurate value of a2 = (91.28 ± 0.2)a.u. 
[stI . I. Blochs group has studied coUisionally driven spin 
dynamics of ®^Rb in Mott insulator regime in an opti- 
cal lattice [1^, and measured very precisely scattering 
length differences. From these measurements we obtain: 
ao = 87.77 ± 0.4a.u., and 04 = 97.23 ± 0.2a.u. That 
implies that 4a4 — 802 = 115. 08a. u. Assuming that 
one may modify spin independent scattering by 10%, 
we get (804 + 4a2)/7 98.88 ± 9.4. We see that 04 
may vary roughly as 04 = 97.28 ± 7.9a.u., along the line 
4a4 — 8a2 = 115. 08a. u. 

This estimate has important consequences: there are 
no feasible OFR, that could reach the regime 04 < ao. 
This is the reason, why reaching the regime of Mott state 
with 2 atoms per site, and total spin 4 (i.e., |2,0,4)) is 
hardly possible with OFRs. 



MPS AND PEPS: A QUANTUM INFORMATION 
APPROACH TO STRONGLY CORRELATED 
SYSTEMS 

The density matrix renormalisation group 
(DMRG)[5i, O is a variational method that has 
had an enormous success in describing ground states 
of some strongly interacting ID systems with rather 
modest computational effort. The underlying philosophy 
of all DMRG oriented algorithms is that many body 
systems can be treated almost " exactly" if one is able to 
truncate the full Hilbcrt space by removing the degrees 
of freedom that arc not involved neither in the ground 



state, nor in the dynamical evolution of the system. The 
difficulty and glory of the method relies on how reliable 



61, m 



the truncation is done. Very recently [55|, |6C 
Quantum Information Theory has provided a new 
perspective on the following questions: (i) how to 
perform an efficient truncation of the Hilbert space, (ii) 
which quantum systems can be efficiently simulated, (iii) 
how to simulate dynamical Hamiltonian and dissipativc 
evolutions of strongly correlated systems, (iv) how and 
when DMRG-orientcd methods can be implemented to 
investigate ground states of 2D and 8D systems; (v) how 
classical concepts like correlation length, which diverge 
on the critical points is linked to the entanglement |64| . 
etc. In perhaps in the simplest version QI approach 
reduces to variational methods based on matrix product 
states, or more general projected entangled pair states In 
general, this Quantum Information approach, apart from 
being extremely simple to implement, is very efficient 
for strongly correlated systems and id already shedding 
a new light in our understanding of many body physics. 
Among recent successes of QI approach are: efficient 
codes for periodic boundary conditions (HBj . simulations 
of finite T and dissipativc systems [63| , renormalization 
algorithms for Quantum-Many Body Systems in two 
and higher dimensions [s^l, understanding of the role 
of entanglement in quantum phase transitions [H^ . 
efficient evaluation of partition functions of frustrated 
and inhomogeneous spin systems [g^, and spectra of 



excited states [22|, studies of quantum impurity modesl 



67| 



[66|, simulation of critical 
quantum lattice systems in ID 
of Laughlin wave functions [b^ 



and infinite-size [68[ 
MPS representations 
simulating quantum 
adiabatic aproach to hard NP-problems [t^, MPS based 
image compressions fn\ . just to name few. 
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